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Abstract 

This chapter deals with charged polymers (polyelectrolytes) in solution and at 
surfaces. The behavior of polyelectrolytes is markedly different from that of neutral 
polymers. In bulk solutions, i.e. disregarding the surface effect, there are two unique 
features to charged polymers: first, due to the presence of long-ranged electrostatic 
repulsion between charged monomers, the polymer conformations are much more 
extended, giving rise to a very small overlap concentration and high solution vis- 
cosity. Second, the presence of a large number of counter-ions increases the osmotic 
pressure of polyelectrolyte solutions, making such polymers water soluble as is of 
great importance to many applications. At surfaces, the same interplay between 
monomer-monomer repulsion and counter-ion degrees of freedom leads to a num- 
ber of special properties. In particular, the adsorption behavior depends on both 
the concentration of polymers and added salt in the bulk. We first describe the 
adsorption behavior of single polyelectrolyte molecules, and discuss the necessary 
conditions to obtain an adsorbed layer and characterize its width. Depending on 
the stiffness of the polyelectrolyte, the layer can be flat and compressed or coiled 
and extended. We then proceed and discuss the adsorption of polyelectrolytes from 
semi-dilute solutions. Mean-field theory profiles of polyelectrolyte adsorption are 
calculated as function of surface charge density (or surface potential), the amount of 
salt in the system and the charge fraction on the chains. The phenomenon of charge 
inversion is reviewed and its relevance to the formation of multilayers is explained. 
The review ends with a short overview of the behavior of grafted polyelectrolytes. 



*to appear in Encyclopedia of Electrochemistry, Eds. M. Urbakh and E. Giladi, Vol. I, Wiley- VCH, 
Weinheim, 2002 
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1 Introduction 



Polyelectrolytes (PE) are charged macromolecules which are extensively studied not only 
because of their numerous industrial applications, but also from a pure scientific inter- 
est [QJ-IQ. The most important property of PE's is their water solubility giving rise to a 
wide range of non-toxic, environmentally friendly and cheap formulations. On the theo- 
retical side, PE's combine the field of statistical mechanics of charged systems with the 
field of polymer science and offer quite a number of surprises and challenges. 

The polymers considered in this review are taken as linear and long polymer chains, 
containing a certain fraction of electrically charged monomers. Chemically, this can be 
achieved, for example, by substituting neutral monomers with acidic ones. Upon con- 
tact with water, the acidic groups dissociate into positively charged protons that bind 
immediately to water molecules, and negatively charged monomers. Although this pro- 
cess effectively charges the polymer molecules, the counter-ions make the PE solution 
electro-neutral on macroscopic length scales. 

The counter-ions are attracted to the charged polymers via long-ranged Coulomb 
interactions, but this physical association typically only leads to rather loosely bound 
counter-ion clouds around the PE chains. Because PE's are present in a background of a 
polarizable and diffusive counter-ion cloud, there is a strong influence of the counter-ion 
distribution on the PE structure, as will be discussed at length in this review. Counter- 
ions contribute significantly towards bulk properties, such as the osmotic pressure, and 
their translational entropy is responsible for the generally good water solubility of charged 
polymers. In addition, the statistics of PE chain conformation is governed by intra-chain 
Coulombic repulsion between charged monomers, resulting in a more extended and swollen 
conformation of PE's as compared to neutral polymers. 

All these factors combined are of great importance when considering PE adsorption 
to charged surfaces. We distinguish between physical adsorption, where chain monomers 
are attracted to surfaces via electrostatic or non-electrostatic interactions, and chemical 
adsorption, where a part of the PE (usually the chain end) is chemically bound (grafted) 
to the surface. In all cases, the long-ranged repulsion of the dense layer of adsorbed PE's 
and the entropy associated with the counter-ion distribution are important factors in the 
theoretical description. 

2 Neutral Polymers in Solution 

Before reviewing the behavior of charged polymers, let us describe some of the important 
ideas underlying the behavior of neutral polymer chains in solution. 

2.1 Flexible Chain Statistics 

The chains considered in this review are either flexible or semi-fiexible. The statistical 
thermodynamics of flexible chains is well developed and the theoretical concepts can 
be applied with a considerable degree of confidence Long and flexible chains 

have a large number of conformations, a fact that plays a crucial role in determining 
their behavior in solution. When flexible chains adsorb on surfaces they form a diffusive 
adsorption layer extending away from the surface into the solution. This is in contrast to 
semi-flexible or rigid chains, which can form dense and compact adsorption layers. 
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The main parameters used to describe a flexible polymer chain are the polymerization 
index A^, which counts the number of repeat units or effective monomers along the chain, 
and the Kuhn length a, being the size of one effective monomer or the distance between 
two neighboring effective monomers. The effective monomer size ranges from a few A 
for synthetic polymers to a few nanometers for biopolymers 0. The effective monomer 
size a is not to be confused with the actual size h of one chemical monomer; in general, 
a > b due to chain stiffening effects as will be explained in detail later on. In contrast 
to other molecules or particles, a polymer chain contains not only translational and rota- 
tional degrees of freedom, but also a vast number of conformational degrees of freedom. 
For typical polymers, different conformations are produced by torsional rotations of the 
polymer backbone bonds. 

A simple description of flexible chain conformations is achieved with the freely-jointed 
chain model where a polymer consisting of + 1 monomers is represented by A^ bonds 
defined by bond vectors with j = 1, . . . A^. Each bond vector has a fixed length \rj\ = a 
corresponding to the Kuhn length, but otherwise is allowed to rotate freely. For the 
freely-jointed chain model, the monomer size b equals the effective monomer size a, b = a. 
Fixing one of the chain ends at the origin, the position of the {k + l)-th monomer is given 
by the vectorial sum 

R, = ^r,. (2.1) 

Because two arbitrary bond vectors are uncorrelated in this simple model, the thermal 
average over the scalar product of two different bond vectors vanishes, (r^ ■ Vk) = for 
j 7^ k, while the mean squared bond vector length is simply given by (r|) = a^. It follows 
then that the mean squared end-to-end radius is proportional to the number of monomers 

(R^) = Na'. (2.2) 

The same result is obtained for the mean quadratic displacement of a freely diffusing 
particle and eludes to the same underlying physical principle, namely the statistics of 
Markov processes. 

In Figure la we show a snapshot of a Monte-Carlo simulation of a freely-jointed chain 
consisting of 100 monomers each being represented by a sphere of diameter b (being equal 
here to a, the effective monomer size). The bar represents a length of 106, which according 
to Eq. ( p.2| ) is the average distance between the chain ends. Indeed, the end-to-end radius 
gives a good idea of the typical chain radius. 

The freely-jointed chain model describes ideal Gaussian chains and does not account 
for interactions between monomers which are not necessarily close neighbors along the 
backbone. Including these interactions will give a different scaling behavior for long 
polymer chains. The end-to-end radius, R = \J {R%), can be written more generally for 
A^ > 1 as 

R ~ aN". (2.3) 

For an ideal polymer chain (no interactions between monomers), the above result implies 
u = 1/2. This result holds only for polymers where the attraction between monomers 
(as compared with the monomer- solvent interaction) cancels the steric repulsion (which 
is due to the fact that the monomers cannot penetrate each other). This situation can be 
achieved in special solvent conditions called "theta" solvents. 
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More generally, polymers in solution can experience three types of solvent conditions, 
with theta solvent condition being intermediates between "good" and "bad" solvent con- 
ditions. The solvent quality depends mainly on the specific chemistry determining the 
interaction between the solvent molecules and monomers. It also can be changed by 
varying the temperature. 

The solvent is called "good" when the monomer-solvent interaction is more favorable 
than the monomer-monomer one. Single polymer chains in good solvents have "swollen" 
spatial configurations, refiecting the effective repulsion between monomers. For good sol- 
vents, the steric repulsion dominates and the polymer coil takes a more swollen structure, 
characterized by an exponent i/ ~ 3/5 0. This spatial size of a polymer coil is much 
smaller than the extended contour length L = aN but larger than the size of an ideal 
chain aN^^^. The reason for this peculiar behavior is entropy combined with the favor- 
able interaction between monomers and solvent molecules in good solvents. As we will see 
below, it is the conformational freedom of polymer coils that leads to salient differences 
between polymer and simple liquid adsorption. 

In the opposite case of "bad" (sometimes called "poor" ) solvent conditions, the effec- 
tive interaction between monomers is attractive, leading to collapse of the chains and to 
their precipitation from the solution (phase separation between the polymer and the sol- 
vent). It is clear that in this case, the polymer size, like any space filling object embedded 
in three dimensional space, scales as ~ i?^, yielding u = 1/3. 

2.2 Semi-Flexible Chain Statistics 

Beside neglecting monomer-monomer interaction, the freely-jointed chain model does not 
take into account the chain elasticity, which plays an important role for some polymers, 
and leads to more rigid structures. This stiffness can be conveniently characterized by 
the persistence length i^, defined as the length over which the tangent vectors at different 
locations on the chain are correlated. In other words, the persistence length gives an 
estimate for the typical radius of curvature, while taking into account thermal fluctuations. 
For synthetic polymers with trans-cis conformational freedom of the chain backbone, the 
stiffness is due to flxed bond angles and hindered rotations around individual back-bone 
bonds. This effect is even more pronounced for polymers with bulky side chains, such as 
poly-DADMAC, because of steric constraints, and the persistence length is of the order 
of a few nanometers. 

Biopolymers with a more complex structure on the molecular level tend to be stiffer 
than simple synthetic polymers. Some typical persistence lengths encountered in these 
systems are £q ^ 5 mm for tubulin, ^ 20 /im for actin, and io ~ 50 nm for double- 
stranded DNA. Because some of these biopolymers are charged, we will discuss at length 
the dependence of the persistence length on the electrostatic conditions. In some cases 
the main contribution to the persistence length comes from the repulsion between charged 
monomers. 

To describe the bending rigidity of neutral polymers, it is easier to use a continuum 
model 0, where one neglects the discrete nature of monomers. The bending energy 
(rescaled by the thermal energy, ksT) of a stiff or semi-flexible polymer of contour length 
L is given by 




(2.4) 
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where d^r(s)/ds^ is the local curvature of the polymer. We assume here that the polymer 
segments are non-extendable, i.e., the tangent vectors are always normalized, |dr(s)/ds| = 
1. Clearly, this continuum description will only be good if the persistence length is larger 
than the monomer size. The mean-squared end-to-end radius of a semi-flexible chain is 
known and reads 

= 2ioL + 2il (e-^/^« - l) , (2.5) 

where the persistence length is io and the total contour length of a chain is L. Two 
limiting behaviors are obtained for R from Eq. ( |2.5| ): for long chains L ^ io, the chains 
behave as flexible ones, R^ ~ 2ioL] while for rather short chains, L <^ io, the chains 
behave as rigid-rods, R c:^ L. Comparison with the scaling of the freely-jointed chain 
model, Eq. ( |2.2| ), shows that a semi-fiexible chain can, for L ^ Iq, be described by a 
freely-jointed chain model with an effective Kuhn length of 

a = 2io (2.6) 

and an effective number of segments or monomers 

In this case the Kuhn length takes into account the chain stiffness and is independent 
from the monomer length. This monomer size is denoted by b whenever there is need to 
distinguish between the monomer size b and the persistence length (or Kuhn length 
a). In Figure 1 we show snapshots taken from a Monte-Carlo simulation of a semi-fiexible 
chain consisting of 100 polymer beads of diameter b. The persistence length is varied from 
io = 2b (Figure lb), over io = 10b (Figure Ic), to io = 1006 (Figure Id). Comparison 
with the freely-jointed chain model (having no persistence length) is given in Figure la 
(a = b, io = 0). It is seen that as the persistence length is increased, the chain structure 
becomes more expanded. The theoretical prediction for the average end-to-end radius R, 
Eq. (p.5|), is shown as the black bar on the figure and gives a good estimate on typical 
sizes of semi-fiexible polymers. 



3 Properties of Polyelectrolytes in Solution 

For polyelectrolytes, electrostatic interactions provide the driving force for their salient 
features and have to be included in any theoretical description. The reduced electrostatic 
interaction between two point-like charges can be written as ziZ2v{r) where 

v(r) = , ^ (3.1) 

is the Coulomb interaction between two elementary charges, zi and Z2 are the valences 
(or the reduced charges in units of the elementary charge e), and e is the medium di- 
electric constant. Throughout this review, all energies are given in units of the thermal 
energy ksT. The interaction depends only on the distance r between the charges. The 
total electrostatic energy of a given distribution of charges is obtained from adding up all 
pairwise interactions between charges according to Eq. (|3.1| ). In principle, the equilibrium 



5 



behavior of an ensemble of charged particles {e.g., a salt solution) follows from the parti- 
tion function, i.e., the weighted sum over all different microscopic configurations, which 
— via the Boltzmann factor — depends on the electrostatic energy of each configuration. 
In practice, however, this route is very complicated for several reasons: 



i) The Coulomb interaction, Eq. (3J), is long-ranged and couples many charged par- 
ticles. Electrostatic problems are typically many-body problems, even for low densities. 

ii) Charged objects in most cases are dissolved in water. Like any material, water is 
polarizable and reacts to the presence of a charge with polarization charges. In addition, 
and this is by far a more important effect, water molecules carry a permanent dipole 
moment that partially orients in the vicinity of charged objects. Note that for water, 
e ~ 80, so that electrostatic interactions and self energies are much weaker in water 
than in air (where e ^ 1) or some other low-dielectric solvents. Still, the electrostatic 
interactions are especially important in polar solvents because in these solvents charges 
dissociate more easily than in unpolar solvents. 

iii) In biological systems and most industrial applications, the aqueous solution con- 
tains mobile salt ions. Salt ions of opposite charge are drawn to the charged object 
and form a loosely bound counter-ion cloud around it. They effectively reduce or screen 
the charge of the object. The effective (screened) electrostatic interaction between two 
charges zic and Z2e in the presence of salt ions and a polarizable solvent can be written 
as ziZ2V-D}i{r) , with the Debye-Hiickel (DH) potential v-onir) given on the linear- response- 
level by 



VBRir) = -e-'^^ (3.2) 



The Bjerrum length £b is defined as 



(3.3) 



and denotes the distance at which the Coulombic interaction between two unit charges in 
a dielectric medium is equal to thermal energy {ksT). It is a measure of the distance be- 
low which the Coulomb energy is strong enough to compete with the thermal fluctuations; 
in water at room temperatures, one finds is ~ 0.7 nm. The exponential decay is charac- 
terized by the so-called screening length which is related to the salt concentration 
Csait by 

= 8nz^iBCs.it (3.4) 

where z denotes the valency of the screening ions. At physiological conditions the salt 
concentration is Cgait ~ 0.1 M and for monovalent ions {z = 1) this leads to k,~^ ~ Inm. 

The so-called Debye-Hiickel interaction, Eq. (pl2|), embodies correlation effects due to 
the long-ranged Coulomb interactions in a salt solution using linear-response theory. [] 
The DH approximation forms a convenient starting point for treating screening effects. 



"'^The number of ions which are correlated in a salt solution with concentration Csait is of the order of 
n ~ K^"^ Csait, where one employs the screening length as the scale over which ions are correlated. 
Using the definition — STrz^^sCgait, one obtains n ~ (^^^sCg^it) The average distance between 

— 1/3 

ions is roughly r ~ c^^^l . The typical electrostatic interaction between two ions in the solution thus is 

U ^ z^is/r ~ ^^^-BCg^it ^'^'i obtain U ^ n~'^/'^ . Using these scaling arguments one obtains that either 
i) many ions are weakly coupled together (i.e., n 3> 1 and U <C 1), or ii) a few ions interact strongly 
with each other ( n ~ [/ ~ 1). In the first case, and in the absence of external fields, the approximations 



leading to the Debye-Hiickel approximation, Eq. (3.2), are valid 
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since (owing to its linear character) the superposition principle is valid and the electro- 
static free energy is given by a sum over the two-body potential Eg. (p.2D. However, we 
will at various points in this review also discuss how to go beyond the DH approximation, 
for example in the form of the non-linear Poisson-Boltzmann theory (see Section ^ or a 
box-model for the counterion distribution (see Section 



3.1 Isolated Polyelectrolyte Chains 

We discuss now the scaling behavior of a single semi-flexible PE in the bulk, including 
chain stiffness and electrostatic repulsion between monomers. For charged polymers, the 
effective persistence length is increased due to electrostatic repulsion between monomers. 
This effect modifies considerably not only the PE behavior in solution but also their 
adsorption characteristics. 

The scaling analysis is a simple extension of previous calculations for flexible (Gaus- 
sian) PE's [|10|-||12[- "^^^ semi-fiexible polymer chain is characterized by a bare persistence 
length io and a linear charge density r. Using the monomer length b and the fraction of 
charged monomers / as parameters, the linear charge density can be expressed as r = f /b. 
Note that in the limit where the persistence length is small and comparable to a monomer 
size, only a single length scale remains, £o — — ^- Many interesting effects, however, are 
obtained in the general case treating the persistence length io and the monomer size b as 
two independent parameters. In the regime where the electrostatic energy is weak, and 
for long enough contour length L, L ^ io, a polymer coil will be formed with a radius R 
unperturbed by the electrostatic repulsion between monomers. According to Eq. (|2.5| ) we 
get ~ 2^0-^- To estimate when the electrostatic interaction will be sufficiently strong 
to swell the polymer coil we recall that the electrostatic energy (rescaled by the thermal 
energy ksT) of a homogeneously charged sphere of total charge Q and radius R is 

1^=, = ^. (3.5) 

5R ^ ' 

The exact charge distribution inside the sphere only changes the prefactor of order unity 
and is not important for the scaling arguments. For a polymer of length L and line charge 
density r the total charge is Q = tL. The electrostatic energy of a (roughly spherical) 
polymer coil is then 

W^ei ^ iBr'L'/%'^'. (3.6) 

The polymer length at which the electrostatic self energy is of order ksT i.e., Wei — 1, 
follows as 

Lei~£o(^BV')~'^' (3.7) 

and defines the electrostatic blob size or electrostatic polymer length. We expect a locally 
crumpled polymer configuration if Lgi > £o, i-e., if 



T^^isio < 1, (3.8) 

because the electrostatic repulsion between two segments of length io is smaller than the 
thermal energy and is not sufficient to align the two segments. This is in accord with more 



detailed calculations by Joanny and Barrat [11]. A recent general Gaussian variational 
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calculation confirms this scaling result and in addition yields logarithmic corrections |T^ . 
Conversely, for 



t^Ib^o > 1, (3.9) 

electrostatic chain-chain repulsion is already relevant on length scales comparable to the 
persistence length. The chain is expected to have a conformation characterized by an 
effective persistence length 4ff, larger than the bare persistence length Eq, i.e., one expects 
4fr > io- 

This effect is clearly seen in Figure 2, where we show snapshots of a Monte-Carlo 
simulation of a charged chain of 100 monomers of size b each and bare persistence length 
io/b = 1 and several values of and r. The number of monomers in an electrostatic 
blob can be written according to Eq. (|3.7| ) as Lei/io = (t^£_b4)^^''^ and yields for Figure 
2a) = 0.25, b) = 0.63, c) Ljio = 1-6, and d) = 4. Accordingly, in d) 

the electrostatic blobs consist of four monomers, and the weakly charged chain crumples 
at small length scales. A typical linear charge density reached with synthetic PE's such 
as Polystyrenesulfonate (PSS) is one charge per two carbon bonds (or, equivalently, one 
charge per monomer), and it corresponds to r 4nm~^. Since for these highly flexible 
synthetic PE's the bare persistence length is of the order of the monomer size, ~ 6, 
the typical charge parameter for a fully charged PE therefore is roughly r'^is^o ~ 3 
and is intermediate between Figures 2a) and 2b). Smaller linear charge densities can 
always be obtained by replacing some of the charged monomers on the polymer backbone 
with neutral ones, in which case the crumpling observed in Figure 2d) becomes relevant. 
Larger bare persistence lengths can be achieved with biopolymers or synthetic PE's with 
a conjugated carbon backbone. 

The question now arises as to what are the typical chain conformations at much larger 
length scales. Clearly, they will be influenced by the repulsions. Indeed, in the persistent 
regime, obtained for Ty/is^o > !> the polymer remains locally stiff even for contour lengths 
larger than the bare persistence length and the effective persistence length is given by 

4ff -4 + W- (3.10) 

The electrostatic persistence length, first derived by Odijk and independently by Skolnick 



and Fixman, reads 13, 14 
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W = ^ (3.11) 

and is calculated from the electrostatic energy of a slightly bent polymer using the lin- 
earized Debye-Hiickel approximation, Eq. (p.2|) . It is valid only for polymer conformation 
which do not deviate too much from the rod-like reference state. The electrostatic per- 
sistence length gives a sizable contribution to the effective persistence length only for 
^OSF > 4- This is equivalent to the condition 



-^/i^o>io'i. (3.12) 



The persistent regime is obtained for parameters satisfying both conditions ( p.9|) and 
(|]T|). Another regime called the Gaussian regime is obtained in the opposite limit of 



The electrostatic persistence length is visualized in Figure 3, where we present snap- 
shots of a Monte-Carlo simulation of a charged chain consisting of 100 monomers of size 
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b. The bare persistence length was fixed at £o = b, and the charge-interaction parame- 
ter was chosen to be t'^C-b^o = 2, close to the typical charge density in experiments on 
fully charged synthetic PE's. The snapshots correspond to varying screening length of a) 
= V2b, leading to an electrostatic contribution to the persistence length of £osf = b, 
b) = \/TSb, or £osF = 96, and c) = V2006, equivalent to £osf = 1006. According 
to the simple scaling principle, Eq. ( |3.1U| ), the effective persistence length in the snap- 
shots. Figure 3a-c, should be similar to the bare persistence length in Figure Ib-d, and 
indeed, the chain structures in 3c) and Id) are very similar. Figure 3a and lb) are 
clearly different, although the effective persistence length should be quite similar, mostly 
because of self-avoidance effects which are present in charged chains and which will be 
discussed in detail in Section p.l.2| . 

For the case where the polymer crumples on length scales larger than the bare bending 
rigidity, i.e., for L^i > io or t^/IbIq < 1, the electrostatic repulsion between polymer 
segments is not strong enough to prevent crumpling on length scales comparable to io, 
but can give rise to a chain stiffening on larger length scales, as explained by Khokhlov 
and Khachaturian and confirmed by Gaussian variational methods [|12[. Figure 4 



schematically shows the PE structure in this regime, where the chain on small scales 
consists of Gaussian blobs of chain length Lei within which electrostatic interactions are 
not important. On larger length scales electrostatic repulsion leads to a chain stiffening, 
so that the PE forms a linear array of electrostatic blobs. To quantify this effect, one 
defines an effective line charge density of a linear array of electrostatic blobs with blob 
size Rei ~ V^oLeh 

~ - ri^^" (3.13) 



Combining Eqs. ( p.l3 ) and ( 3.11 ) gives the effective electrostatic persistence length for a 



string of electrostatic blobs, 

()l/3^-2/3^2/3 



'B ^0 



(3.14) 



This electrostatic stiffening is only relevant for the so-called Gaussian-persistent regime 
valid for £kk > -Rei, or equivalently 



"V/^ > (4«:)'/'. (3.15) 



When this inequality is inverted the Gaussian persistence regime crosses over to the 
Gaussian one. 

The crossover boundaries (|3.9|), ( |3.12| ), ( |3.15| ) between the various scaling regimes are 



summarized in Figure 5. We obtain three distinct regimes. In the persistent regime, for 
ta/^bA) > and ta/^b^o > 1) the polymer takes on a rod-like structure with an effective 
persistence length larger than the bare persistence length and given by the OSF expression, 
Eq. ( |3.11| ). In the Gaussian-persistent regime, for ta/Tb^o < 1 and t\/Tb^q > (ioK,)^^'^, the 



polymer consists of a linear array of Gaussian electrostatic blobs, as shown in Figure 4, 
with an effective persistence length £kk larger than the electrostatic blob size and given 
by Eq. (|3.14|) . Finally, in the Gaussian regime, for ta/Jb^q < i^of^Y^'^ and t^/Ib^ < ^o^; 
the electrostatic repulsion does not lead to stiffening effects at any length scale. 

The persistence length £kk was also obtained from Monte-Carlo simulations with pa- 
rameters shown in Figure 2d), where chain crumpling at small length scales and chain 
stiffening at large length scales occur simultaneously |1^, [1^, |1^, |T^. However, extremely 



9 



long chains are needed in order to obtain reliable results for the persistence length, since 
the stiffening occurs only at intermediate length scales and therefore fitting of the tangent- 
tangent correlation function is nontrivial. Nevertheless, simulations point to a different 
scaling than in Eg. ( |3.14] ), with a dependence on the screening length closer to a lin- 
ear one, in qualitative agreement with experimental results^. The situation is com- 
plicated by the fact that more recent theories for the single polyelectrolyte chain give 
different results, some confirming the simple scaling results described in Eqs.( p.ll]) and 



(CT).P. 0, some confirming Eq.(pil) while criticizing Eq.(P7[^P, |2T|, This 



issue is not resolved and under intense current investigation. For multivalent counterions 
fluctuation effects can even give rise to a polyelectrolyte collapse purely due to electro- 



static interactions [23, 24, 25|, which is accompanied by a negative contribution to the 



effective persistence length p6 



3.1.1 Manning Condensation 

A peculiar phenomenon occurs for highly charged PE's and is known as the Manning 



condensation of counter- ions |l27| , Pq] , For a rigid PE represented by an infinitely long and 



straight cylinder with a linear charge density larger than 

isTZ = 1 (3.16) 

where z is the counter-ion valency, it was shown that counter-ions condense on the oppo- 
sitely charged cylinder even in the limit of infinite solvent dilution. Real polymers have a 
finite length, and are neither completely straight nor in the infinite dilution limit [^, |3^ . 



Still, Manning condensation has an experimental significance for polymer solutions be- 
cause thermodynamic quantities, such as counter-ion activities and osmotic coeffi- 
cients [^, show a pronounced signature of Manning condensation. Locally, polymer 



segments can be considered as straight over length scales comparable to the persistence 
length. The Manning condition Eq. (|3.16|) usually denotes a region where the binding of 
counter-ions to charged chain sections begins to deplete the solution from free counter- 
ions. Within the scaling diagram of Figure 5, the Manning threshold (denoted by a 
vertical broken line) is reached typically for charge densities larger than the one needed 
to straighten out the chain. This holds for monovalent ions provided £o > ^s? as is almost 
always the case. The Manning condensation of counter-ions will not have a profound 
influence on the local chain structure since the chain is rather straight already due to 
monomer-monomer repulsion. A more complete description of various scaling regimes 
related to Manning condensation, chain collapse and chain swelling has recently been 



given [33 . 



3.1.2 Self- Avoidance and Polyelectrolyte Chain Conformations 

Let us now consider how the self-avoidance of PE chains comes into play, concentrating 
on the persistent regime defined by ta/^^^o > 1- The end-to-end radius i? of a strongly 
charged PE chain shows three distinct scaling ranges. For a chain length L smaller than 
the effective persistence length igg, which according to Eq. ( |3.10|) is the sum of the bare and 



electrostatic persistence lengths, R grows linearly with the length, R ^ L. Self-avoidance 
plays no role in this case, because the chain is too short to fold back on itself. 

For much longer chains, L ^ 4fr, we envision a single polymer coil as a solution of 
separate polymer pieces of length fgfr, and treat their interactions using a virial expansion. 
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The second virial coefficient f 2 of a rod of length ies and diameter d scales as f 2 ~ i'^gd jT^ . 
The chain connectivity is taken into account by adding the entropic chain elasticity as a 
separate term. The standard Flory theory for a semi-flexible chain is based on writing 
the free energy JF as a sum of two terms 



^ d3 / ^/^' 



2 



+ (3.17) 



i?3 



where the ffist term is the entropic elastic energy associated with swelling a polymer chain 
to a radius R and the second term is the second-virial repulsive energy proportional to 
the coefficient V2 and the segment density squared. It is integrated over the volume R^. 
The optimal radius R is calculated by minimizing this free energy and gives the swollen 
radius 

Rr^ iv2/ iesY^'L" (3.18) 

with p = 3/5. This swollen radius is only realized above a minimal chain length L > 
Lsw ~ ^Iff/'^2 ~ ^es/^"^- For elongated segments with i^s 3> d, or, equivalently, for a 
highly charged PE, we obtain an intermediate range of chain lengths < L < L^w for 
which the chain is predicted to be Gaussian and the chain radius scales as 

R^i%^L^/\ (3.19) 

For charged chains, the effective rod diameter d is given in low salt concentrations by 
the screening length, i.e., d ~ k,~^ plus logarithmic corrections. The condition to have 
a Gaussian scaling range, Eq. ( p.l9| ), thus becomes 4fr 3> Within the persistent 

and the Gaussian-persistent scaling regimes depicted in Figure 5 the effective persistence 
length is dominated by the electrostatic contribution and given by Eqs. (|3.11|) and ( p.14 ) 



respectively, which in turn are always larger than the screening length It follows that 
a Gaussian scaling range, Eq. ( p.l9|) , always exists below the asymptotic swollen scaling 
range, Eq. ( |3.18| ). This situation is depicted in Figure 4 for the Gaussian-persistent scaling 



regime, where the chain shows two distinct Gaussian scaling ranges at the small and large 
length scales. This multihierarchical scaling structure is only one of the many problems 
one faces when trying to understand the behavior of PE chains, be it experimentally, 
theoretically, or by simulations. 

A different situation occurs when the polymer backbone is under bad-solvent condi- 
tions, in which case an intricate interplay between electrostatic chain swelling and short- 



range collapse occurs|^^. Quite recently, this interplay was theoretically shown to lead 
to a Rayleigh instability in the form of a necklace structure consisting of compact beads 
connected by thin strings P5|, [57| , Small-angle X-ray scattering on solvophobic 



polyelectrolytes in a series of polar organic solvents of various solvent quality could qual- 
itatively conffim these theoretical predictions []39[| . 



3.2 Dilute Polyelectrolyte Solutions 

It is natural to generalize the discussion of single chain behavior to that of many PE 
chains in dilute solvent concentrations. The dilute regime is defined by Cm < cj^, where 
Cm denotes the monomer concentration (per unit volume) and is the concentration 
where individual chains start to overlap. Clearly, the overlap concentration is reached 
when the average bulk monomer concentration exceeds the monomer concentration inside 
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a polymer coil. To estimate the overlap concentration c^, we simply note that the average 
monomer concentration inside a coil with radius R ~ bN" is given by 

C^^^ N'-'^h-\ (3.20) 

For ideal chains with z/ = 1/2 the overlap concentration scales as ~ N"^^"^ and thus 
decreases slowly as the polymerization index N increases. For rigid polymers with v = 1 
the overlap concentration scales as c*^ ~ iV"^ and decreases strongly as increases. This 
means that the dilute regime for semi-flexible PE chains corresponds to extremely low 
monomer concentrations. For example taking a monomer size h = 0.254 nm and a poly- 
merization index of = 10^, the overlap concentration becomes ~ 6 x 10" ''nm-^ ^ 
10~3 mM, which is a very small concentration. 

The osmotic pressure in the dilute regime in the limit Cm — >^ is given by 



ksT z N 

and consists of the ideal pressure of non-interacting counter-ions (first term) and polymer 
coils (second term). Note that since the second term scales as A^~^, it is quite small for 
large N and can be neglected. Hence, the main contribution to the osmotic pressure comes 
from the counter-ion entropy. This entropic term explains also why charged polymers can 
be dissolved in water even when their backbone is quite hydrophobic. Precipitation of the 
PE chains will also mean that the counter-ions are confined within the precipitate. The 
entropy loss associated with this confinement is too large and keeps the polymers dispersed 
in solution. In contrast, for neutral polymers there are no counter- ions in solution. Only 
the second term in the osmotic pressure exists and contributes to the low osmotic pressure 
of these polymer solutions. In addition, this can explain the trend towards precipitation 
even for very small attractive interactions between neutral polymers. 

3.3 Semi-Dilute Polyelectrolyte Solution 

In the semi-dilute concentration regime, > c^, different polymer coils are strongly over- 
lapping, but the polymer solution is still far from being concentrated. This means that 
the volume fraction of the monomers in solution is much smaller than unity, b^Cm <C 1. In 
this concentration range, the statistics of counter-ions and polymer fluctuations are inti- 
mately connected. One example where this feature is particularly prominent is furnished 
by neutron and X-ray scattering from semi-dilute PE solutions [^-^^. The structure 



factor S{q) shows a pronounced peak, which results from a competition between the con- 
nectivity of polymer chains and the electrostatic repulsion between charged monomers, as 
will be discussed below. An important length scale, schematically indicated in Figure 6, 
is the mesh-size or correlation length C,, which measures the length below which entan- 
glement effects between different chains are unimportant. The mesh size can be viewed 
as the polymer (blob) scale below which single-chain statistics are valid. A semi-dilute 
solution can be thought of being composed of roughly a close-packed array of polymer 
blobs of size ^. 

The starting point for the present discussion is the screened interaction between 
two charges immersed in a semi-dilute PE solution containing, charged polymers, their 
counter-ions and, possibly, additional salt ions. Screening in this case is produced not 
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only by the ions, but also by the charged chain segments which can be easily polarized 
and shield any free charges. 

Using the random-phase approximation (RPA), the effective Debye-Hiickel (DH) in- 
teraction can be written in Fourier space as p6l HI] 



/ N 1 + V2CmSo{q) 

recalling that Cm is the average density of monomers in solution and / is the fraction of 
charged monomers on each of the PE chains. The second virial coefficient of monomer- 
monomer interactions is V2 and the single-chain form factor (discussed below) is denoted 
by 5*0 (g). In the case where no chains are present, = 0, the RPA expression reduces to 
'^rpa(q') = VDuiq), the Fourier-transform of the Debye-Hiickel potential of Eq. (^.2]), given 
by 

^DH(g) = 4t^2- (3-23) 

q'' + K'^ 

As before, is the DH screening length, which is due to all mobile ions. We can write 
^2 = Kga^it + AnzlBfCm^ whcre Kga^n- = Sh z'^ i BCsait dcscribcs the screening due to added salt 
of concentration Cgait? and the second term describes the screening due to the counter- 
ions of the PE monomers. Within the same RPA approximation the monomer-monomer 
structure factor S{q) of a polymer solution with monomer density Cm is given by ^ 

S^\q) = Pv^n{q) + S^\q)/c^ + V2. (3.24) 

The structure factor (or scattering function) only depends on the form factor of an iso- 
lated, non- interacting polymer chain, 5*0 (g), the second virial coefficient, f2, the fraction / 
of charged monomers, and the interaction between monomers, which in the present case is 
taken to be the Debye-Hiickel potential fDH('?)- The structure factor of a non- interacting 
semi-flexible polymer is characterized, in addition to the monomer length 6, by its persis- 
tence length £eff- In general, this form factor is a complicated function which cannot be 



written down in closed form However, one can separate between three different 

ranges of wavenumbers g, and within each range the form factor shows a rather simple 
scaling behavior, namely 



N'^ for q^<Q/NM, 



cS 



So\q)-{ q^hU/Q for 6/iV64ff < < 36/7rW (3.25) 
gfo/vr for 36/7r2£2^ < q^. 

For small wavenumbers the polymer acts like a point scatterer, while in the intermediate 
wavenumber regime the polymer behaves like a flexible, Gaussian polymer, and for the 
largest wavenumbers the polymer can be viewed as a stiff rod. 

One of the most interesting features of semi-dilute PE solutions is the fact that the 
structure factor S{q) shows a pronounced peak. For weakly charged PE's, the peak 
position scales as g ~ with the monomer density in agreement with the above 



random-phase approximation (RPA) results for charged polymers P6| , We now discuss 
the scaling of the characteristic scattering peak within the present formalism. The position 
of the peak follows from the inverse structure factor, Eq. ( ^.24| ) , via dS~^ {q)/dq = 0, which 
leads to the equation 

,V.L.+4..../c„^(|g^)'^ (3.26) 
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In principle, there are two distinct scaling behaviors possible for the peak, depending 
on whether the chain form factor of Eq. ( |3.25|) exhibits flexible-like or rigid-like scaling. 



Concentrating now on the flexible case, i.e. the intermediate g-range in Eq. (|3.25|) , the 
peak of S{q) scales as 

in agreement with experimental results. A peak is only obtained if the left-hand side of 
Eq. ( |3.26| ) is dominated by the g-dependent part, i.e., if {q*Y > /tg^^n + AnziBfCm- 



In Figure 7a we show density-normalized scattering curves for a PE solution charac- 
terized by the persistence length £es = 1 nm (taken to be constant and thus independent 
of the monomer concentration), with monomer length b = 0.38 nm (as appropriate for 
poly-DADMAC) and charge fraction / = 0.5 and with no added salt. As the monomer 
density decreases (top to bottom in the figure), the peak moves to smaller wavenumbers 
and sharpens, in agreement with previous implementations of the RPA. In Figure 7b 
we show the same data in a different representation. Here we clearly demonstrate that 
the large-g region already is dominated by the 1/q behavior of the single-chain structure 
factor, So{q). Since neutron scattering data easily extend to wavenumbers as high as 
q ~ 5nm~^, the stiff-rod like behavior in the high g-limit, exhibited on such a plot, will 
be important in interpreting and fitting experimental data even at lower g-values. 

In a semi-dilute solution there are three different, and in principle, independent length 
scales: The mesh size ^, the screening length k~^, and the persistence length ies- In the 
absence of added salt, the screening length scales as 

~ (ziBfcm)-'^^ . (3.28) 

Assuming that the persistence length is larger or of the same order of magnitude as the 
mesh size, as is depicted in Figure 6, the polymer chains can be thought of straight 
segments between different crossing links. Denoting the number of monomers inside a 
correlation blob as g, this means that ^ ^ bg. The average monomer concentration scales 
as Cm ^ g/^^ , from which we conclude that 

e~(&c™)-'/'. (3.29) 

Finally, the persistence length within a semi-dilute PE solution can be calculated by con- 
sidering the electrostatic energy cost for slightly bending a charged rod. In PE solutions, 
it is important to include in addition to the screening by salt ions also the screening 
due to charged chain segments. This can be calculated by using the RPA interaction, 
Eq. ( p.22| ). Since the screening due to polymer chains is scale dependent and increases 
for large separations, a g-dependent instability is encountered and leads to a persistence 
length 

^OSF ~ ibcmr'^' (3.30) 

where the 'sd' superscript stands for 'semi-dilute'. This result is a generalization of the 
OSF result for a single chain and applies to semi-dilute solutions. Comparing the three 
lengths, we see that 

e ~ ^OSF - (3.31) 
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Since the prefactor \J C-bI /b for synthetic fully charged polymers is roughly of order unity, 
one finds that for salt-free semi-dilute PE solutions, all three length-scales scale in the 
same way with c^, namely as ~ c~^/^, as is known also from experiments pO| , p| , |5T| and 



previous theoretical calculations p2|, p^ ]. In simulations of many polyelectrolyte chains, the 
reduction of the chain size due to screening by PE chains was clearly seen|]5^, |55|, p^, |57 . 



4 Adsorption of a Single Polyelectrolyte Chain 

After reviewing bulk properties of PE solutions we address the complete adsorption di- 
agram of a single semi- flexible polyelectrolyte on an oppositely charged substrate. In 
contrast to the adsorption of neutral polymers, the resulting phase diagram shows a large 
region where the adsorbed polymer is flattened out on the substrate and creates a dense 
adsorption layer. 

Experimentally, the adsorption of charged polymers on charged or neutral substrates 
has been characterized as a function of the polymer charge, chemical composition of the 



substrate, pH and ionic strength of the solution [58, B9[, as well as the substrate charge 



density Repeated adsorption of anionic and cationic PE's can lead to well char- 

acterized multilayers on planar [^-|^ and spherical substrates |70|. Theoretically, 



the adsorption of PE's on charged surfaces poses a much more complicated problem than 
the corresponding adsorption of neutral polymers. The adsorption process results from 
a subtle balance between electrostatic repulsion between charged monomers, leading to 
chain stiffening, and electrostatic attraction between the substrate and the polymer chain. 
The adsorption problem has been treated theoretically employing the uniform expansion 
method [|7l| and various continuous mean-field theories [j72|-|]75[. In all these works, the 



polymer density is taken to be constant in directions parallel to the surface. 

The adsorption of a single semi-flexible and charged chain on an oppositely charged 
plane can be treated as a generalization of the adsorption of flexible polymers . 



A PE characterized by the linear charge density r, is subject to an electrostatic potential 
created by a, the homogeneous surface charge density (per unit area). Because this 
potential is attractive for an oppositely charged substrate, we consider it as the driving 
force for the adsorption. More complex interactions are neglected. They are due to the 
dielectric discontinuity at the substrate surface and to the impenetrability of the substrate 
for salt ions. 

Within the linearized DH theory the electrostatic potential of a homogeneously charged 
plane is 

W(a^) = 47r£Bfx«:-ie-«^ (4.1) 

Assuming that the polymer is adsorbed over a layer of width 5 smaller than the screening 
length the electrostatic attraction force per monomer unit length can be written as 

/att = -^nlBcrr (4.2) 

We neglect non-linear effects due to counter-ion condensation. They are described by the 
Gouy-Chapman theory for counter-ion distribution close to a charged surface. Although 



^An ion in solution has a repulsive interaction from the surface when the solution dielectric constant 
is higher than that of the substrate. This effect can lead to desorption for highly charged PE chains. On 
the contrary, when the substrate is a metal there is a possibility to induce PE adsorption on non-charged 
substrates or on substrates bearing charges of the same sign as the PE. See ref. [Tq for more details. 
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these effects are clearly important, it is difficult to include them systematically, and we 
remain at the linearized Debye-Hiickel level. 

Because of the confinement in the adsorbed layer, the polymer feels an entropic repul- 
sion. If the layer thickness S is much smaller than the effective persistence length of the 
polymer, i^s, as depicted in Figure 8a, a new length scale, the so-called deflection length 
A, enters the description of the polymer statistics. The deflection length A measures the 
average distance between two contact points of the polymer chain with the substrate. As 
shown by Odijk, the deflection length scales as A ~ 6'^^^il(^ and is larger than the layer 



thickness S but smaller than the persistence length ies The entropic repulsion follows 
in a simple manner from the deflection length by assuming that the polymer loses roughly 
a free energy of one ksT per deflection length. 

On the other hand, if ^ > ies, as shown in Figure 8b, the polymer forms a random coil 
with many loops within the adsorbed layer. For a contour length smaller than L ~ 
the polymer obeys Gaussian statistics and decorrelates into blobs with an entropic cost 
of one ksT per blob. The entropic repulsion force per monomer unit length is thus ITS 



f r^l ^''^Xs^' for '^«4fr ^4 3^ 

^'"^ I ies6-' for 5»4fr ^ ' 

where we neglected a logarithmic correction factor which is not important for the scaling 
arguments. As shown in the preceding section, the effective persistence length 4fr depends 
on the screening length and the line charge density; in essence, one has to keep in mind 
that £eG is larger than io for a wide range of parameters because of electrostatic stiffening 
effects, g 

The equilibrium layer thickness follows from equating the attractive and repulsive 
forces, Eqs. ( [4.2|) and ( [4.3|) . For rather stiff polymers and small layer thickness, 6 < 
< ies, we obtain 

6 {isarill,')''^' (4.4) 

For a layer thickness corresponding to the screening length, 6 ~ k^^, scaling arguments 
predict a rather abrupt desorption transition . This is in accord with previous Transfer- 



Matrix calculations for a semi-flexible polymer bound by short-ranged (square-well) po- 
tentials [^, |80|, Setting S ~ in Eq. (|4^ ), we obtain an expression for the 



adsorption threshold (for nics > 1) 

^5/3 

a* ^ (4.5) 

For a > a* the polymer is adsorbed and localized over a layer with a width smaller than 
the screening length (and with the condition ies > also satisfying 6 < 4ff)- As a 
is decreased, the polymer abruptly desorbs at the threshold a = a* . In the Gaussian 
regime, the effective persistence length ies is given by the bare persistence length £o and 

•^The situation is complicated by the fact that the electrostatic contribution to the persistence length 
is scale dependent and decreases as the chain is bent at length scales smaller than the screening length. 



This leads to modifications of the entropic confinement force, Eq. { \i.3\ j , if the defiection length is smaller 
than the screening length. As can be checked explicitly, all results reported here are not changed by these 
modifications. 
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the desorption threshold is obtained by replacing 4ff by io in Eq. (^.51), i.e., 



5/3 

a* (4.6) 



In the persistent regime, we have ics ~ ^osf with Iqsf given by Eq. (|3.11|) . The adsorption 
threshold follows from Eq. (|4.5| ) as 

Finally, in the Gaussian-persistent regime, we have an effective line charge density from 
Eq. ( p.l3| ) and a modified persistence length, Eq. ( |3.14| ). For the adsorption threshold we 
obtain from Eq. ( [4.5|) 

^7/3/<5/9 

Let us now consider the opposite limit, £efr < Q If the layer thickness is larger 
than the persistence length but smaller than the screening length, £efr < S < k~^, the 
prediction for S obtained from balancing Eqs. ([4. 21) and (|4.3|) becomes 



(. 



f \ 1/3 
f-cff \ 



(4.9) 



From this expression we see that 6 has the same size as the screening length k^^ for 

a* ~ (4.10) 

This in fact denotes the location of a continuous adsorption transition at which the layer 
grows to infinity. The scaling results for the adsorption behavior of a flexible polymer, 
Eqs. ( |4.9| )-( ^mi| ), are in agreement with previous results [j7T|. 

In Figure 9 we show the desorption transitions and the line at which the adsorbed 
layer crosses over from being flat, 6 < 4fr, to being crumpled or coiled, 6 > ics- The 
underlying PE behavior in the bulk, as shown in Figure 5, is denoted by broken lines. We 
obtain two different phase diagrams, depending on the value of the parameter 

S = ail^^e]^^ (4.11) 

For strongly charged surfaces, S > 1, we obtain the phase diagram as in Figure 9a, and 
for weakly charged surfaces, S < 1, as in Figure 9b. We see that strongly charged PE's, 
obeying rv%?B > 1, always adsorb in fiat layers. The scaling of the desorption transitions 
is in general agreement with recent computer simulations of charged PE's P^j. 
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From Eq. (4.4) we see that the layer thickness S is of the same order as £cs for isO'Ti'^g ~ 1, at which 



point the condition 5 <^ i^s used in deriving Eq. (4.4) breaks down 
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4.1 Adsorption on Curved Substrates 



Adsorption of polyelectrolytes on curved substrates is of importance because PE's are 
widely used to stabilize colloidal suspensions ||83[ and to fabricate hollow polymeric 
shells |7D[. When the curvature of the small colloidal particles is large enough, it 
can lead to a much more pronounced effect for PE adsorption as compared with neutral- 
polymer. This is mainly due to the fact that the electrostatic energy of the adsorbed PE 
layer depends sensitively on curvature p4l RR p6|. Bending a charged polymer around a 



small sphere costs a large amount of electrostatic energy, which will disfavor adsorption 
of long, strongly charged PE at too low salt concentration. 

In Figure 10 we show the adsorption phase diagram of a single stiff PE of finite length 
which interacts with an oppositely charged sphere of charge Z (in units of e). The specific 
parameters were chosen as appropriate for the complexation of DNA (a negatively charged, 
relatively stiff biopolymer) with positively charged histone proteins, corresponding to a 
DNA length of L = 50 nm, a chain persistence length of £o = 30 nm, and a sphere radius 
of i?p = 5nm. The phase diagram was obtained by minimization of the total energy 
including bending energy of the DNA, electrostatic attraction between the sphere and 
the DNA, and electrostatic repulsion between the DNA segments with respect to the 
chain configuration ||8^. Configurational fluctuations away from this ground state are 
unimportant for such stiff polymers. 

We show in Figure 10 the main transition between an unwrapped state, at low sphere 
charge Z, and the wrapped state, at large sphere charge Z. It is seen that at values of 
the sphere charge between Z = 10 and Z = 130 the wrapping only occurs at intermediate 



values of the inverse screening length k 



.1/2 
-'salt 



. At low salt concentrations, (lower left cor- 



ner in the phase diagram), the self-repulsion between DNA segments prevents wrapping, 
while at large salt concentrations, (lower right corner in the diagram), the electrostatic 
attraction is not strong enough to overcome the mechanical bending energy of the DNA 
molecule. These results are in good agreement with experiments on DNA/histone com- 
plexes |^8|. Interestingly, the optimal salt concentration, where a minimal sphere charge 
is needed to wrap the DNA, occurs at physiological salt concentration, for f» Inm. 
For colloidal particles of larger size and for flexible synthetic polymers, configurational 
fluctuations become important. They have been treated using a mean-field description in 



terms of the average monomer density profile around the sphere [89, 90 



5 Adsorption from Semi-Dilute Solutions 

So far we have been reviewing the behavior of single PE chains close to a charged wall 
(or surface). This will be now extended to include adsorption of PE from bulk (semi- 
dilute) solutions having a bulk concentration c^. As before the chains are assumed to 
have a fraction / of charged monomers, each carrying a charge e, resulting in a linear 
charge density r = //6 on the chain. The solution can also contain salt (small ions) of 
concentration Cgait which is directly related to the Debye-Hiickel screening length, k,~^. 
For clarity purposes, the salt is assumed to be monovalent {z = 1) throughout Section 5. 

We will consider adsorption only onto a single flat and charged surface. Clearly the 
most important quantity is the profile of the polymer concentration Cm{x) as function of 
X, the distance from the wall. Another useful quantity is the polymer surface excess per 
unit area, defined as 
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Cm X 



(5.1) 



Related to the surface excess V is the amount of charges (in units of e) carried by the 
adsorbing PE chains, fV. In some cases the polymer carries a higher charge (per unit 
area) than the charged surface itself, fV > a, and the surface charge is overcompensated 
by the PE as we will see later. This does not violate charge neutrality in the system 
because of the presence of the counter-ions in solution. 

In many experiments, the total amount of polymer surface excess P is measured as a 
function of the bulk polymer concentration, pH and/or ionic strength of the bulk solution 



g and 



)1[-|9^. For reviews see, e.g., refs. |5^, |100|| . More recently, spectroscopy 
ellipsometry |^ have been used to measure the width of the adsorbed PE layer. Other 
techniques such as neutron scattering can be employed to measure the entire profile Cm{x) 
of the adsorbed layer |p.01| , p.02 



In spite of the difficulties to treat theoretically PE's in solution because of the deli- 
cate interplay between the chain connectivity and the long range nature of electrostatic 
interactions [|I|, |103| , |104|| , several simple approaches treating adsorption exist. One 
approach is a discrete multi-Stern layer model ||105|| - ||T09| , where the system is placed on 
a lattice whose sites can be occupied by a monomer, a solvent molecule or a small ion. 
The electrostatic potential is determined self-consistent ly (mean-field approximation) to- 
gether with concentration profiles of the polymer and small ions. In another approach, 
the electrostatic potential and the PE concentration are treated as continuous functions 
ITII , |73| , |74| , |110| , |111| , |112| , |1 13|| . These quantities are obtained from two coupled differen- 
tial equations derived from the total free energy of the system. We will review the main 
results of the latter approach, presenting numerical solutions as well as scaling arguments 
of the mean field profiles. 



5.1 Mean-Field Theory and Its Profile Equations 

The charge density on the polymer chains is assumed to be continuous and uniformly 
distributed along the chains. Further treatments of the polymer charge distribution (an- 
nealed and quenched models) can be found in refs. |^4|, |1 1 1|| . Within mean-field approx- 
imation, the free energy of the system can be expressed in terms of the local electrostatic 
potential V'(r), the local monomer concentration Cm(r) and the local concentration of pos- 
itive and negative ions c^(r). The mean-field approximation means that the infiuence of 
the charged surface and the inter-chain interactions can be expressed in term of an exter- 
nal potential which will determine the local concentration of the monomers, Cm,(r)- This 
external potential depends both on the electrostatic potential and on the excluded vol- 
ume interactions between the monomers and the solvent molecules. The excess free energy 
with respect to the bulk can then be calculated using another important approximation, 
the ground state dominance. This approximation is used often for neutral polymers |^ 
and is valid for very long polymer chains, ^ 1. It is then convenient to introduce the 
polymer order parameter 0(r) where Cm(r) = |0(r)P and to express the adsorption free 
energy JF in terms of and (and in units of ksT) |7^, |7^, |110| , |111| , |112| 



J dr {Fpoi(r) + Fions(r) + Fei(r)} . (5.2) 
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The polymer contribution is 



Fpoi(r) = y|V0r + ^t;2(0^-0^)-M0'-0') , (5.3) 

where the first term is the polymer elastic energy. Throughout this section we restrict 
ourselves to flexible chains and treat the Kuhn length a and the effective monomer size h 
as the same parameter. The second term is the excluded volume contribution where the 
second virial coefficient V2 is of order a?. The last term couples the system to a polymer 
reservoir via a chemical potential yUp, and (p^ = \Jd^ is related to the bulk monomer 
concentration, c^. 

The entropic contribution of the small (monovalent) ions is 

i^ions(r) = 1^ d -d - Csalt In Csalt + Csalt " Ai*(c* - Cgalt) , (5.4) 

i=± 

where c*(r) and /i* are, respectively, the local concentration and the chemical potential of 
the i = ± ions, while Csait is the bulk concentration of salt. 
Finally, the electrostatic contributions (per ksT) are 

Fei(r) = [/e02^ + ec+^-ec-^-|^|VV^|2| AbT . (5.5) 

The first three terms are the electrostatic energies of the monomers, the positive ions 
and the negative ions, respectively, / is the fractional charge carried by one monomer. 
The last term is the self energy of the electric field where e is the dielectric constant of 
the solution. Note that the electrostatic contribution, Eq. (pT5|), is equivalent to the well 
known result: (e/STrfc^T) / dr IVV'P plus surface terms. This can be seen by substituting 
the Poisson-Boltzmann equation (as obtained below) into Eq. ( |5.5| ) and then integrating 
by parts. 

Minimization of the free energy, Eqs. ( p.2[ )-( p75D with respect to c^, and ip yields 
a Boltzmann distribution for the density of the small ions, c^(r) = Csait exp(=Fe-?/'//cBT), 
and two coupled differential equations for (p and ip: 

V'^ir) = — Cait sinh(e^/fc^T) - — (f<j>' - /^^e^^/'^-^) , (5.6) 

^V20(r) = t;2(0=^-0g0) + /0e7/>/fcBT , (5.7) 





Equation (|5.6| ) is a generalized Poisson-Boltzmann equation including the free ions as 



well as the charged polymers. The first term represents the salt contribution and the 
second term is due to the charged monomers and their counter-ions. Equation (|5.7| ) is a 
generalization of the self-consistent field equation of neutral polymers In the bulk, 
the above equations are satisfied by setting -^/^ — and (pb- 



5.2 Numerical Profiles: Constant t/jg 

When the surface is taken as ideal, i.e.,, fiat and homogeneous, the physical quanti- 
ties depend only on the distance x from the surface. The surface imposes boundary 
conditions on the polymer order parameter (j){x) and electrostatic potential tp{x). In ther- 
modynamic equilibrium all charge carriers in solution should exactly balance the surface 
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charges (charge neutrahty). The Poisson-Bohzmann equation ( |5.6| ), the self-consistent 
field equation ( ^.71) and the boundary conditions uniquely determine the polymer concen- 
tration profile and the electrostatic potential. In most cases, these two coupled non-linear 
equations can only be solved numerically. 

We present now numerical profiles obtained for surfaces with a constant potential ips'- 

^L=o = ^- (5-8) 

The boundary conditions for depend on the nature of the short range non-electrostatic 
interactions of the monomers and the surface. For simplicity, we take a non-adsorbing 
surface and require that the monomer concentration will vanish there: 

0L=o = . (5.9) 



We note that the boundary conditions chosen in Eqs. (|5.8|)-(|5l^) model the particular 
situation of electrostatic attraction in competition with a short range (steric) repulsion of 
non-electrostatic origin. Possible variations of these boundary conditions include surfaces 
with a constant surface charge (discussed below) and surfaces with a non-electrostatic 
short range attractive (or repulsive) interaction with the polymer [|7^, |115|| . Far from 



the surface {x oo) both ip and reach their bulk values and their derivatives vanish: 

ii'\x^oo = and (p'lx^oo = 0. 

The numerical solutions of the mean-field equations (|5.6|), ( p.7[ ) together with the 
boundary conditions discussed above are presented in Figure 11, for several different 
physical parameters. The polymer is positively charged and is attracted to the non- 
adsorbing surface held at a constant negative potential. The aqueous solution contains 
a small amount of monovalent salt (cgait = 0.1 mM). The reduced concentration profile 
Cm{x)/(f)l is plotted as a function of the distance from the surface x. Different curves 
correspond to different values of the reduced surface potential i/s = eips/ksT, the charge 
fraction / and the effective monomer size a. Although the spatial variation of the profiles 
differs in detail, they all have a single peak which can be characterized by its height and 
width. This observation serves as a motivation to using scaling arguments. 



5.3 Scaling Results 

The numerical profiles of the previous section indicate that it may be possible to obtain 
simple analytical results for the PE adsorption by assuming that the adsorption is char- 
acterized by one dominant length scale D. Hence, we write the polymer order parameter 
profile in the form 

0(a;) = ^h{x/D) , (5.10) 

where h{x) is a dimensionless function normalized to unity at its maximum and cm sets 
the scale of polymer adsorption, such that 0(-D) = ^/cM■ The free energy can now be 
expressed in terms of D and cm while the exact form of h{x) affects only the numerical 
prefactors. 

In principle, the adsorption length D depends also on the ionic strength through k,~^. 
As discussed below the scaling assumption, Eq. ( |5.10| ), is only valid as long as and 
D are not of the same order of magnitude. Otherwise, h should be a function of both kx 
and x/D. We concentrate now on two limiting regimes where Eq. ( |5.10| ) can be justified: 
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(i) the low-salt regime D <^ and (ii) the high-salt regime D ^ /t^^. We first discuss 
the case of constant surface potential which can be directly compared with the numerical 
profiles. Then we note the differences with the constant surface charge boundary condition 
where the interesting phenomenon of charge overcompensation is discussed in detail. 



5.3.1 Low— Salt Regime D <^ k and = constant 

In the low-salt regime the effect of the small ions can be neglected and the free energy 
(per unit surface area), Eqs. (p.2|)-( p3D , is approximated by (see also refs. [|73], p,12|| ) 



F ^ ^CM - flVslcMD + ATTlBfcl,D' + ^V2CI,D . (5.11) 

In the above equation and in what follows we neglect additional prefactors of order unity 
in front of the various terms which arise from inserting the scaling profile Eq. ( |5.10|) into 
the free energy. The first term of Eq. (|5.11|) is the elastic energy characterizing the 



response of the polymer to concentration inhomogeneities. The second term accounts 
for the electrostatic attraction of the polymers to the charged surface. The third term 
represents the Coulomb repulsion between adsorbed monomers. The last term represents 
the excluded volume repulsion between adsorbed monomers, where we assume that the 
monomer concentration near the surface is much larger than the bulk concentration Cm ^ 

In the low-salt regime and for highly charged PE's the electrostatic interactions are 
much stronger than the excluded volume ones. Neglecting the latter interactions and 
minimizing the free energy with respect to D and Cm gives: 

and 



|2 



CM - ^-\^s\\ (5.13) 

recalling that ys = eips/ksT. As discussed above, these expressions are valid as long as 
(i) D ^ K,^^ and (ii) the excluded volume term in Eq. ( [5.11| ) is negligible. Condition (i) 



translates into Csait f\ys\/{S7ilBCi^)- For \ys\ ~ 1, a = 5A and Ib = 7A this limits the 
salt concentration to Csa.it/ f 0.4 M. Condition (ii) on the magnitude of the excluded 
volume term can be shown to be equivalent to / ^ V2\ys\/J'B0''^- These requirements are 
consistent with the data presented in Figure 11. 

We recall that the profiles presented in Figure 11 were obtained from the numerical 
solution of Eqs. ( ^.6|) and (|5.71 ), including the effect of small ions and excluded volume. 
The scaling relations are verified by plotting in Figure 12 the same sets of data as in 
Figure 11 using rescaled variables as defined in Eqs. (|5.12|) , (|5.13|) . Namely, the rescaled 



electrostatic potential ip{x)/ips and polymer concentration Cm{x)/cM ~ Cm{x)a /\ys\ are 
plotted as functions of the rescaled distance x/D ~ xf^^'^\ys\^^'^ /a. The different curves 
roughly collapse on the same curve. 

In many experiments the total amount of adsorbed polymer per unit area F is measured 
as function of the physical characteristics of the system such as the charge fraction /, the 
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pH of the solution or the salt concentration Csait [^Il'fE^- This quantity can be easily 
obtained from our scaling expressions yielding 

rco \y |3/2 \Jj 13/2 

r = [c^{x) - <Pl]dx ^ DcM j^j^^ ~ . (5.14) 

The adsorbed amount r(/) in the low-salt regime is plotted in the inset of Figure 13a. 
As a consequence of Eq. ( ^.14|) , F decreases with increasing charge fraction /. Similar 



behavior was also reported in experiments |^4[. This effect is at first glance quite puzzling 



because as the polymer charge increases, the chains are subject to a stronger attraction 
to the surface. On the other hand, the monomer-monomer repulsion is stronger and 
indeed, in this regime, the monomer-monomer Coulomb repulsion scales as (/ca/)^, and 
dominates over the adsorption energy which scales as /cm- 

5.3.2 High-Salt Regime D ^ k^^ and = constant 

Let us now consider the opposite case of a high ionic strength solution. Here D is much 
larger than and the electrostatic interactions are short ranged with a cut-off k~^. 
The free energy of the adsorbing PE layer (per unit surface area) then reads: 

F ~ -^cm - f\ys\cMK-^ + A'kIbPk-^I.D + \v2cl,D . (5.15) 

The electrostatic cut-off enters in two places. In the second term only the first layer of 
width interacts electrostatically with the surface. In the third term each charged layer 
situated at point x interacts only with layers at x' for which |a; — a;'| < kT^ . This term 
can be also viewed as an additional electrostatic excluded volume with fgi ~ Isif /t^Y- 
Minimization of the free energy gives 

,2 1/2 



and 



yielding 



f\ys\ fws\ 



f\ys\Csli'^ /IV^sl ^-1/2 , . 



where a is a numerical constant of order unity which depends on the profile details. 

The adsorption behavior is depicted in Figures 13 and 14. Our results are in agreement 
with numerical solutions of discrete lattice models (the multi-Stern layer theory) [|58|, 
H, |105|, |10|, |03, [10|, |10|. In Figure 13, F is plotted as function of / (Figure 13a) and 



the pH (Figure 13b) for different salt concentrations. The behavior as seen on Figure 13b 
represents annealed PE's where the nominal charge fraction is given by the pH of the 
solution through the expression 

IQpH-pKo 
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where pKg = — log^g and Kq is the apparent dissociation constant. We note that this 
relation is only strictly valid for infinitely dilute monomers and that distinct deviations 
from it are observed for polyelectrolytes because of the electrostatic repulsion between 
neighboring dissociating sites ||114(| . Still, the results in Figure 13b capture the main 



qualitative trends of pH-dependent PE adsorption. Another interesting observation which 

1/2 

can be deduced from Eq. ( |5.18| ) is that T is only a function of fc^^^il ■ Indeed, as can be 
seen in Figure 13, Csait only affects the position of the peak and not its height. 

The effect of salt concentration is shown in Figure 14, where F is plotted as function of 
the salt concentration Csait for two charge fractions / = 0.01 and 0.25. The curves on the 
right hand side of the graph are calculated from the high-salt expression for F, Eq. ( FTsI) . 



The horizontal lines on the left hand side of the graph indicate the low-salt values of F, 
Eq. (|5.14|) . The dashed lines in the intermediate salt regime serve only as guides to the 
eye since our scaling approach is not valid when D and k,~^ are of the same order. 

Emphasis should be drawn to the distinction between weakly and strongly charged 
PE's. For weak PE's, the adsorbed amount F is a monotonously decreasing function of 
the salt concentration Csait in the whole range of salt concentrations. The reason being 
that the monomer-monomer Coulomb repulsion, proportional to is weaker than the 
monomer-surface interaction, which is linear in /. 

For strongly charged PE's, on the other hand, the balance between these two electro- 
static terms depends on the amount of salt. At low salt concentrations, the dominant 
interaction is the monomer-surface Coulomb repulsion. Consequently, addition of salt 
screens this interaction and increases the adsorbed amount. When the salt concentration 
is high enough, this Coulomb repulsion is screened out and the effect of salt is to weaken 
the surface attraction. At this point the adsorbed amount starts to decrease. As a result, 
the behavior over the whole concentration range is non-monotonic with a maximum at 
some optimal value c^^^^ as seen in Figure 14. 

From this analysis and from Figures 13,14 and Eq. ( ^.18 ), it is now natural to divide 



the high-salt regime into two sub regimes according to the PE charge. At low charge 
fractions (sub-regime HS I), / ^ /* = {csa.itV2Y^'^, the excluded volume term dominates 
the denominator of Eq. ( |5.18|) and 

r ~ /lV^s|c-f , (5.20) 

whereas at high / (sub-regime HS II), / ^ /*, the monomer-monomer electrostatic 
repulsion dominates and F decreases with / and increases with Csait^ 

r ~ cll^Mr'- (5.21) 

The various regimes with their crossover lines are shown schematically in Figure 15. 
Keeping the charge fraction / constant while changing the amount of salt corresponds to 
a vertical scan through the diagram. For weak PE's this cut goes through the left hand 
side of the diagram starting from the low-salt regime and, upon addition of salt, into the 
HS I regime. Such a path describes the monotonous behavior inferred from Figure 14 for 
the weak PE (/ = 0.01). For strong PE the cut goes through the right hand side of the 
diagram, starting from the low-salt regime, passing through the HS II regime and ending 
in the HS I. The passage through the HS II regime is responsible to the non-monotonous 
behavior inferred from Figure 14 for the strong PE (/ = 0.25). 

Similarly, Figures 13a, 13b correspond to horizontal scans through the top half of the 
diagram. As long as the system is in the HS I regime, the adsorbed amount increases when 
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the polymer charge fraction increase. As the polymer charge further increases, the system 
enters the HS II regime and the adsorbed amount decreases. Thus, the non-monotonous 
behavior of Figure 13. We finally note that the single-chain desorption transition for a 
flexible chain which was discussed in Section 4, is also valid for adsorption from solutions 
and will lead to a desorption transition at very high salt concentrations in Fig. 15. 



5.4 Overcompensation of Surface Charges: Constant a 

We turn now to a different electrostatic boundary condition of constant surface charge 
density and look at the interesting phenomenon of charge compensation by the PE chains 
in relation to experiments for PE adsorption on fiat surfaces, as well as on charged colloidal 



particles |T0|. What was observed in experiments is that PE's adsorbing on an 

oppositely charged surface can overcompensate the original surface charge. Because the 
PE's create a thin layer close to the surface, they can act as an effective absorbing surface 
to a second layer of PE's having an opposite charge compared to the first layer. Repeating 
the adsorption of alternating positively and negatively charged PE's, it is possible to create 
a multilayer structure of PE's at the surface. Although many experiments and potential 
applications for PE multilayers exist, the theory of PE overcompensation is only starting 



to be developed 0, |TT2|, [TT3HTT5|, [TT6| . 



The scaling laws presented for constant %jjs can be used also for the case of constant 
surface charge. A surface held at a constant potential ips will induce a surface charge 
density a. The two quantities are related by: dip/dx = inae/e at x = 0. We will now 
consider separately the two limits: low salt D ^ and high salt D ^ n^^. 



5.4.1 Low Salt Limit: D < k'^ 

Assuming that there is only one length scale characterizing the potential behavior in the 
vicinity of the surface, as demonstrated in Figure 12a, we find that the surface potential 
ips and the surface charge a are related by ips ~ creD. In the low salt limit we find from 
Eq.(CT) 

D ~ if^rhr'^' . (5.22) 

Let us define two related concepts via the effective surface charge density defined as 
Act = /r — cr, which is sum of the adsorbed polymer charge density and the charge density 
of the bare substrate. For Aa = the adsorbed polymer charge exactly compensates the 
substrate charge. If Aa is positive the PE overcompensates the substrate charge, more 
polymer adsorbs than is needed to exactly cancel the substrate charge. If Aa is positive 
and reaches the value Aa = a it means that the PE charge is /F = 2a and leads to a 
charge inversion of the substrate charge. The effective surface charge consisting of the 
substrate charge plus the PE layer has a charge density which is exactly opposite to the 
original substrate charge density a. 

Do we obtain overcompensation or charge inversion in the low salt limit within mean- 
field theory? Using scaling arguments this is not clear since we find that Aa ~ /F ~ 
a. Namely each of the two terms in Aa scales linearly with a, and the occurrence of 
overcompensation or charge inversion will depend on numerical prefactors determining 
the relative sign of the two opposing terms. However, if we look on the numerical solution 
for the mean- field electrostatic potential. Figure 12, we see indeed that all plotted profiles 
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have a maximum of ^(x) as function of x. An extremum in means a zero local electric 
field. Or equivalently, using Gauss law, this means that the integrated charge density 
from the wall to this special extremum point (including surface charges) is exactly zero. 
At this point the charges in solution exactly compensate the surface charges. 



5.4.2 High Salt Limit: D > /t^ 

When we include salt in the solution and look at the high salt limit the situation is more 
complex. The only length characterizing the exponential decay of ijj close to the surface 
is the Debye-Hiickel screening length. Hence, using dip/dx\s ~ aeK~^ yields ips ~ cre/t"^ 
and therefore from Eq. (|5.16|) 

Dr^^r. K^f-'a-' (5.23) 

The estimation of the PE layer charge can be obtained by using the expression for D and 
Cm in this high salt limit, Eqs. (|5.16|) - (|5.18|) , yielding 



^ /3a(87r/BCsait/t ^) ^ (3a 

1 + af2Csalt/~^ l + V2/Vel 

where Vei = /^/acgait is the electrostatic contribution to the second virial coefficient V2, 
and a > and f3 > 1 are positive numerical factors. 

We see that Aa = /F — a is a decreasing function of V2- Charge overcompensation 
can occur when V2 is smaller than Ve\ (up to a prefactor of order unity). When V2 can be 
neglected in the vicinity of the surface, or when f 2 = (theta solvents), there is always 
charge overcompensation. Act = (/? — l)a > 0. This is the case of strongly charged PE. 



Similar conclusions have been mentioned in refs. ||75|, |115|| where f 2 was taken as zero but 
the surface has a non-electrostatic short range interaction with the PE. By tuning the 
relative strength of the surface charge density a and the non-electrostatic interaction it is 
also possible to cause a charge overcompensation and even an exact charge inversion in a 
special case. 

Finally, we note that the dependence of the charge parameter Aa on the amount of 
salt, Csait, is different for constant surface charge and constant surface potential cases. 
While for the former Act is non-monotonous and has a maximum (as mentioned above) 
as function of the salt concentration, in the latter case Aa ja is a monotonic decreasing 
function of Csait, Eq. (|5.24]) . This can be explained by the extra powers of Csait in the latter 
case coming from the relation ^/^^ ~ ~ aK~^ . 

Let us remark that in other theories the overcharging is due to lateral correlations 
between adsorbed PE's, which in conjunction with screening by salt ions leads to strongly 



overcharged surfaces [76, 116 . 



5.5 Final Remarks on Adsorption from Semi-Dilute Solutions 

The results presented above (Section 5) for adsorption from solutions have been derived 
using mean-field theory. Hence, lateral fluctuations in the polymer and ionic concentra- 
tions are neglected. In addition, we neglect the delicate influence of the charges on the 
PE persistence length and any deviations from ground state dominance 0. The region 
of validity of the theory is for long and weakly charged polymer chains in contact with a 
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moderately charged surface. The PE solution is placed in contact with a single and ideal 
surface (infinite, flat and homogeneous). The problem reduces then to an effective one 
dimensional problem depending only on the distance from the charged surface. We take 
very simple boundary conditions for the surface assuming that the polymer concentration 
is zero on the surface and keeping the surface in constant potential or constant surface 
charge conditions. 

We find numerical solutions for the PE profile equations in various cases. These 
numerical solutions agree well with simple scaling assumptions describing the adsorption 
of PE's. Scaling expressions for the amount of adsorbed polymer F and the width D of 
the adsorbed layer, as a function of the fractional charge / and the salt concentration 
Csait, are obtained for two cases: constant iIjs and constant a. 

For constant ifjg and in the low-salt regime a /"^^^ dependence of F is found. It 
is supported by our numerical solutions of the profile equations (|5.6|), ( |5.7| ) and is in 



agreement with experiment |^4[. This behavior is due to strong Coulomb repulsion 
between adsorbed monomers in the absence of salt. As / decreases, the adsorbed amount 
increases until the electrostatic attraction becomes weaker than the excluded volume 
repulsion, at which point F starts to decrease rapidly. At high salt concentrations it 
is not possible to neglect the excluded volume interaction of the monomers since the 
electrostatic interactions are screened by the salt. We obtain two limiting behaviors: (i) 
for weakly charged PE's, / ^ /* = (csait'i^2)^''^, the adsorbed amount increases with the 
fractional charge and decreases with the salt concentration, F ~ //i/cWt, due to the 
monomer-surface electrostatic attraction, (ii) For strong PE's, / ^ /*, the adsorbed 
amount decreases with the fractional charge and increases with the salt concentration, 
r ~ v^Csait/ /, due to the dominance of the monomer-monomer electrostatic repulsion. 
Between these two regimes we find that the adsorbed amount reaches a maximum in 
agreement with experiments . 



The scaling arguments are then repeated for constant o boundary conditions. It is 
found that the PE can possibly cause charge overcompensation and even inversion of 
the nominal substrate charge, leading the way to multilayer formation of positively and 
negatively charged PE. The scaling approach can serve as a starting point for further 
investigations. The analytical and approximated expressions are valid only in specific 
limits. Special attention should be directed to the crossover regime where D and are 
of comparable size. 

The problem of charge inversion is not well understood at present. Alternative ap- 
proaches rely on lateral correlations between semi-fiexible adsorbed PE chains which also 



can lead to strong overcompensation of surface charges [76, 116 



6 Polyelectrolyte Brushes 



Charged polymers which are densely end-grafted to a surface are called polyelectrolyte 
or charged brushes. They have been the focus of numerous theoretical |TT^-|T55[ and 
experimental ||126| , |127| , |128|| studies. In addition to the basic interest, charged brushes are 
considered for their applications as efficient means for preventing colloids in polar media 
(such as aqueous solutions) from flocculating and precipitating out of solution ||83|. This 
stabilization arises from steric (entropic) as well as electrostatic (energetic) repulsion. 
A strongly charged brush is able to trap its own counter-ions and generates a layer of 
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locally enhanced salt concentration ||119||. It is thus less sensitive to the salinity of the 



surrounding aqueous medium than a stabilization mechanism based on pure electrostatics 
{i.e., without polymers). 

Neutral brushes have been extensively studied theoretically in the past using scaling 
theories | 129 , 130 1, strong-stretching theories ||131| , 132 , 133 1, self-consistent field theo- 



ries 



134||, and computer simulations 



Little is known from experiments on the 
scaling behavior of PE brushes as compared to uncharged polymer brushes. The thickness 
of the brush layer has been calculated from neutron-scattering experiments on end-grafted 
polymers ||126|| and charged diblock-copolymers at the air- water interface |p.28| . 

Theoretical work on PE brushes was initiated by the works of Miklavic and Marcelja 
117|| and Misra et al. [|118|| . In 1991, Pincus ||119|| and Borisov, Birshtein and Zhulina 



120|| presented scaling theories for charged brushes in the so-called osmotic regime, where 



the brush height results from the balance between the chain elasticity (which tends to 
decrease the brush height) and the repulsive osmotic counter-ion pressure (which tends 
to increase the brush height). In later studies, these works have been generalized to poor 
solvents ||121|| and to the regime where excluded volume effects become important, the 
so-called quasi- neutral or Alexander regime | 124 ]. 



In what follows we assume that the charged brush is characterized by two length 
scales: the average vertical extension of polymer chains from the wall L, and the typical 
extent of the counter-ion cloud, denoted by H. We neglect the presence of additional 
salt, which has been discussed extensively in the original literature, and only consider 
screening effects due to the counter-ions of the charged brush. Two different scenarios 
emerge, as is schematically presented in Figure 16. The counter-ions can either extend 
outside the brush, H ^ L, as shown in 16a), or be confined inside the brush, H ^ L, as 
shown in 16b). As we show now, case b) is indicative of strongly charged brushes, while 
case a) is typical for weakly charged brushes. 

The free energy per unit area (and in units of ksT) contains several contributions. 
We denote the grafting density of PE's by p, the counter-ion valency by z, recalling that 

is the polymerization index of grafted chains and / their charge fraction. The osmotic 
free energy, Fqs, associated with the ideal entropy cost of confining the counter-ions to a 
layer of thickness H is given by 



TP ~ 

-'OS — 



Nfp^^fNfp' 
z \ zH 



(6.1) 



is the second virial contribution to the free energy, arising from steric repulsion 
between the monomers (contributions due to counter ions are neglected). Throughout 
this review, the polymers are assumed to be in a good solvent (positive second virial 
coefficient V2 > 0). The contribution thus reads 



- 2 



(6.2) 



Finally, a direct electrostatic contribution Fgi occurs if the PE brush is not locally electro- 
neutral throughout the system, as for example is depicted in Figure 16a. This energy is 
given by 

^ 27ciB{Nfp)^L-Hf 



H 



(6.3) 
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This situation arises in the hmit of low charge, when the counter-ion density profile extends 
beyond the brush layer, i.e., H > L. 

The last contribution is the stretching energy of the chains which is 

Here, a is the monomer size or Kuhn length of the polymer, implying that we neglect any 
chain stiffness for the brush problem. The different free energy contributions lead, upon 
minimization with respect to the two length scales H and L, to different behaviors. Let 
us first consider the weak charging limit, i.e., the situation where the counter-ions leave 
the brush, H > L. In this case, minimization of F^s + F^i with respect to the counter-ion 
height H leads to 

which is the Gouy-Chapman length for 2;-valent counter-ions at a surface of surface charge 
density a = N f p. Balancing now the polymer stretching energy F^^ and the electrostatic 
energy F^i one obtains the so-called Pincus brush 

Lc^N^pa^lsf ■ (6.6) 

In the limit of if ~ L the PE brush can be considered as neutral and the electrostatic 
energy vanishes. There are two ways of balancing the remaining free energy contributions. 
The first is obtained by comparing the osmotic energy of counter- ion confinement, Fqs, 
with the polymer stretching term, Fgt, leading to the height 



Naf 



1/2 



,1/2 ' (6-7) 



constituting the so-called osmotic brush regime. Finally comparing the second-virial free 
energy, F^j, with the polymer stretching energy, Fst, one obtains 

L Na{v2p/af''^ , (6.8) 

and the PE brush is found to have the same scaling behavior as the neutral brush ||129| , 



130|| . Comparing the heights of all three regimes we arrive at the phase diagram shown 



in Figure 17. The three scaling regimes meet at the characteristic charge fraction 



and the characteristic grafting density 

For large values of the charge fraction / and the grafting density p it has been found numer- 
ically that the brush height does not follow any of the scaling laws discussed here ||137| . 



This has been recently rationalized in terms of another scaling regime, the collapsed 
regime. In this regime one finds that correlation and fluctuation effects, which are ne- 
glected in the discussion in this section, lead to a net attraction between charged monomers 



and counter-ions 1138 
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7 Conclusions 



In this chapter we have reviewed the behavior of charged polymers (polyelectrolytes) in 
solution and at interfaces, concentrating on aspects that are different from the correspond- 
ing behavior of neutral polymers. 

Because polyelectrolytes (PE) tend to be quite stiff due to electrostatic monomer- 
monomer repulsions, their chain statistics is related to that of semi-flexible polymers. 
Neutral and charged semi-flexible polymers are controlled by their bending rigidity, which 
is usually expressed in terms of a persistence length (see Section 2.2). For PE's, the 
electrostatic interaction influences considerably this persistence length. 

In solution, we have considered the scaling behavior of a single PE (Section 3.1). The 
importance of the electrostatic persistence length was stressed. The Manning condensa- 
tion of counter-ions leads to a reduction of the effective linear charge density (Section 
3.1.1). Excluded volume effects are typically less important than for neutral polymers 
(Section 3.1.2). Dilute PE solutions are typically dominated by the behavior of the 
counter-ions. So is the large osmotic pressure of dilute PE solutions due to the entropic 
contribution of the counter-ions (Section 3.2). Semi-dilute PE solutions can be described 
by the random-phase-approximation (RPA), which in particular yields the characteristic 
peak of the structure factor. 

At surfaces, we discussed in detail the adsorption of single PE (Section 4), the adsorp- 
tion from semi-dilute solutions (Section 5), and the behavior of end-grafted PE chains 
(Section 6). We tried to express the PE behavior in terms of a few physical parame- 
ters such as the chain characteristics (persistence length), ionic strength of the solutions 
and surface characteristics. The shape and size of the adsorbing layer is, in many in- 
stances, governed by a delicate balance of competing mechanisms of electrostatic and 
non-electrostatic origin. In some cases it is found that the adsorbing PE layer is flat and 
compressed, while in other cases it is coiled and extended. Yet, in other situations, the 
PE will not adsorb at all and will be depleted from the surface. We also briefly review 
the phenomenon of charge overcompensation and inversion, when the adsorbed PE layer 
effectively inverses the sign of the surface charge leading the way to formation of PE 
multilayers. 

Important topics that we have left out are the dynamics of PE solutions, which is 
reviewed in ref. Q, and the behavior of PE under bad-solvent conditions |36, 37, 38[ . 



In the future we expect that studies of PE in solutions and at surfaces will be directed 
more towards biological systems. We mentioned in this review the complexation of DNA 
and histones (Sec. 4.1). This is only one of many examples of interest where charged 
biopolymers, receptors, proteins and DNA molecules interact with each other or with 
other cellular components. The challenge for future fundamental research will be to try 
to understand the role of electrostatic interactions combined with specific biological (lock- 
key) mechanisms and to infer on biological functionality of such interactions. 
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Legend of Symbols 

• a: Kuhn length (equivalent to monomer size for flexible polymers) 

• b: monomer size in the case of semi-flexible polymers 

• Cm' monomer concentration in dilute solution. 

• Cm{x): monomer profile at distance x from the surface 

• c^: bulk monomer concentration in semi-dilute solutions. 

• c^: overlap concentration of bulk polymer solution 

• Cm' characteristic value of PE adsorption concentration (semi-dilute case) 

• Csait^ salt concentration in the solution 

• c^(a;): profiles of ± ions 

• d: charged rod diameter (or cross-section) 

• D: adsorption layer thickness (semi-dilute case) 

• S: thickness of adsorption layer (single chain case) 

• e: electronic unit charge 

• e: dielectric constant of the medium. £ = 80 for water. 

• /: fraction charge of the chain < / < 1 

• /*: crossover value for / 

• fatt- attractive (electrostatic) adsorption force (per monomer unit length) (single 
chain) 

• /rep: repulsive (entropic) adsorption force (per monomer unit length) (single chain) 

• excess free energy of adsorption per unit area (semi-dilute) 

• Fpoi- polymer contribution to the free energy density (semi-dilute) 

• -^ions^ entropic contribution of small ions (semi-dilute) 

• Fei: electrostatic contribution (semi-dilute) 

• Fog-, osmotic counter- ion contribution to brush free energy 

• Fy2: second- virial contribution to brush free energy 

• Fei: direct electrostatic contribution to brush free energy 

• Fgt: chain stretching contribution to brush free energy 
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• 4>{x) = ^Cm{x): polymer order parameter 

• (f)},: bulk value of polymer order parameter 

• g: number of monomer per blob 

• F: polymer surface excess per unit area 

• h{x): dimensionless PE adsorption profile 

• H: counter-ion cloud vertical extend from the surface (PE brush case) 

• /cfiT: thermal energy 

• Debye-Hiickel screening length 

• Ksait^ salt contribution to k 

• L: contour length of a chain 

• L: polymer vertical extent from the surface (PE brush case) 

• Lei: chain length inside one electrostatic blob 

• Lg^: chain length inside one swollen blob (Sec. 3.1.2) 

• ^B — e^/{£kBT): Bjerrum length 

• £o- bare (mechanical) persistence length 

• ^osp: electrostatic contribution to persistence length (Odijk, Skolnick, and Fixman 
length) 

• ^OSF- generahzation of the OSF to semi-dilute regime. 

• 4ff: effective persistence length 

• ^kr: persistence length of a string of blobs 

• //^: chemical potential of ± ions 

• ^p. chemical potential of polymer 

• N: polymerization index 

• I/: Flory exponent for the polymer size 

• H: osmotic pressure 

• %l){x): electrostatic potential at point x from the surface 

• ■05: surface potential at x = 

• q*: the q wavenumber at the peak of S{q) 

• Q: total charge of particle or polyelectrolyte chain 
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• R: end-to-end polymer chain radius 

• it!ei: radius of one electrostatic blob 

• Rp-. spherical particle radius (Sec. 4.1) 

• Rjv: end-to-end chain vector 

• Tj: position vector of the ith monomer 

• p: grafting density of a PE brush 

• p*: crossover value for p 

• s: < s < L continuous contour chain parameter 

• So{q): form factor of a single chain 

• S{q): structure factor (or scattering function) of a PE solution 

• a: surface charge density (in units of e) at x = 

• a*: minimal threshold value for a for single chain adsorption 

• Act — fr — a: overcharging parameter 

• T, — atj'^t^: dimensionless adsorption parameter 

• T = f /h: linear charge density on the chain 

• V2: 2nd virial coefficient of monomers in solution. ^2 > for good solvents. 

• v{r) = I {kBTer): dimensionless Coulomb interaction between two ions 

• v-£)Yi(yr) = v{r) exp{—Kr): Debye-Hiickel interaction 

• vbb.{q)'- DH interaction in Fourier space 

• ^rpa(9): RPA interaction in Fourier space 

• y{x) = eip{x)/kBT: dimensionless potential profile 

• Us — y(0): rescaled surface potential 

• ^: correlation length (mesh size) of semi-dilute polymer solution 

• z — ±1, ±2, . . . valency of the ions 

• Z: charge of the spherical particle in units of e (Sec. 4.1) 
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Figure Captions 



Figure 1: Snapshots of Monte-Carlo simulations of a neutral and semi-fiexible chain 
consisting of = 100 monomers of diameter b, which defines the unit of length. The 
theoretical end-to-end radius R is indicated by a straight bar. The persistence lengths 
used in the simulations are: a) io = 0, corresponding to a freely-jointed (flexible) chain, 
leading to an end-to-end radius R/b = 10, b) io/b = 2, leading to R/b = 19.8, c) io/b = 10, 
leading to R/b = 42.4, d) 4/6 = 100, leading to R/b = 85.8. 

Figure 2: Snapshots of Monte-Carlo simulations of a polyelectrolyte chain of = 100 
monomers of size b, taken as the unit length. In all simulations the bare persistence 
length is fixed at io/b = 1, and the screening length and the charge interactions are 
tuned such that the electrostatic persistence length is constant and iosp/b = 100. The 
parameters used are: a) K~^/b = ^/EO and t'^Hb^q = 8, b) /b = V200 and r'^is^o = 2, 
c) K-^/b = VSOO and r^is^o = 1/2, and d) n'^b = ^3200 and rHnio = 1/8. Noticeably, 
the weakly charged chains crumple at small length scales and show a tendency to form 
electrostatic blobs. 

Figure 3: Snapshots of Monte-Carlo simulations of a PE chain consisting of = 100 
monomers of size b. In all simulations, the bare persistence length is fixed at io/b = 1, and 
the charge-interaction parameter is chosen to be t^^b^o = 2. The snapshots correspond 
to varying screening length of: a) /b = \/2, leading to an electrostatic contribution 
to the persistence length of iosp/b = 1, b) /b = a/TS, leading to iosp/b = 9, and c) 
K^^/b = v200, leading to iosp/b = 100. According to the simple scaling principle, Eq. 
( ^.1U| ), the effective persistence length in the snapshots a)-c) should be similar to the bare 
persistence length in Figure lb)-d). 

Figure 4: Schematic view of the four scaling ranges in the Gaussian-persistent regime. 
On spatial scales smaller than R^i the chain behavior is Gaussian; on length scales larger 
than i?ei but smaller than £kk the Gaussian blobs are aligned linearly. On length scales 
up to Lsvj the chain is isotropically swollen with an exponent = 1/2, and on even larger 
length scales self-avoidance effects become important and u changes to z/ = 3/5. 

Figure 5: Schematic phase diagram of a single semi-fiexible PE in bulk solution with bare 
persistence length £q and line charge density r, exhibiting various scaling regimes. High 
salt concentration and small r correspond to the Gaussian regime, where the electrostatic 
interactions are irrelevant. In the persistent regime, the polymer persistence length is 
increased, and in the Gaussian-persistent regime the polymer forms a persistent chain 
of Gaussian blobs as indicated in Figure 4. The broken lines indicates the Manning 
condensation, at which counter-ions condense on the polymer and reduce the effective 
polymer line charge density. We use a log-log plot, and the various power-law exponents 
for the crossover boundaries are denoted by numbers. 

Figure 6: Schematic view of the chain structure in the semi-dilute concentration range. 
The mesh size ^ is about equal to the effective polymer persistence length i^s and to the 
screening length (if no salt is added to the system). 



Figure 7: a) RPA prediction for the rescaled structure factor S{q)/cm of a semi-dilute 
PE solution with persistence length iea = Inm, monomer length b = 0.38 nm and charge 
fraction / = 0.5 in the salt-free case. The monomer densities are (from bottom to top) 
Cm =1M, 0.3 M, 10 mM, 3mM, ImM, and 0.3 mM. b) For the same series of values as 
in a) the structure factor is multiplied by the wavenumber q. The semi-flexibility becomes 
more apparent because for large q the curves tend towards a constant. 

Figure 8: a) Schematic picture of the adsorbed polymer layer when the effective persis- 
tence length is larger than the layer thickness, ies > S. The distance between two contacts 
of the polymer with the substrate, the so-called deflection length, scales as A ~ 6'^^^il(^. 
b) Adsorbed layer for the case when the persistence length is smaller than the layer thick- 
ness, ies < In this case the polymer forms a random coil with many loops and a 
description in terms of a flexible polymer model becomes appropriate. 

Figure 9: Adsorption scaling diagram shown on a log-log plot for a) strongly charged 
surfaces, S = ai^^'^i^^'^ > 1 and for b) weakly charged surfaces S < 1. We find a desorbed 
regime, an adsorbed phase where the polymer is flat and dense, and an adsorbed phase 
where the polymer shows loops. It is seen that a fully charged PE is expected to adsorb 
in a flat layer, whereas charge-diluted PE's can form coiled layers with loops and dangling 
ends. The broken lines denote the scaling boundaries of PE chains in the bulk as shown 
in Figure 5. The numbers on the lines indicate the power law exponents of the crossover 
boundaries between the regimes. 

Figure 10: Numerically determined adsorption diagram for a charged semi-flexible poly- 
mer of length L = 50 nm, linear charge density r = 6 nm~^, persistence length = 30 nm, 
interacting with an oppositely charged sphere of radius Rp = 5nm. Shown is the main 
transition from the unwrapped configuration (at the bottom) to the wrapped configuration 
(at the top) as a function of sphere charge Z and inverse screening length k. Wrapping is 
favored at intermediate salt concentrations. The parameters are chosen for the problem 
of DNA-histone complexation. Adapted from ref. 



Figure 11: Adsorption profiles obtained by numerical solutions of Eqs. ( p. 61) , ( ^.7|) for 
several sets of physical parameters in the low-salt limit. The polymer concentration scaled 
by its bulk value = (pl is plotted as a function of the distance from the surface. The 
different curves correspond to: / = 1, a = 5A and i/s = eipg/ksT = —0.5 (solid curve); 
/ = 0.1, a = 5A and i/s = —0.5 (dots); / = 1, a = 5A and i/s = —1.0 (short dashes); 
f = 1, a = lOA and i/s = —0.5 (long dashes); and / = 0.1, a = 5A and i/s = 1.0 (dot-dash 
line). For all cases 0f = lO'^A-^, V2 = 50A^, e = 80, T = 300K and Csait = 0.1 mM. 



Adapted from ref. ||112 



Figure 12: Scaling behavior of PE adsorption in the low-salt regime Eqs. ( p.l2| ), ( p.l3| ). 



(a) The rescaled electrostatic potential il^{x)/\il^s\ as a function of the rescaled distance 
x/D. (b) The rescaled polymer concentration Cm{x)/cM as a function of the same rescaled 
distance. The profiles are taken from Figure 11 (with the same notation). The numerical 
pref actors of a linear h{x/D) profile were used in the calculation of D and cm- Adapted 
from ref. ||112|| . 



Figure 13: Typical adsorbed amount F as a function of (a) the charge fraction / and 
(b) the pH — pKq of the solution for different salt concentrations Eq. ( ^.18| ). The insets 



correspond to the low-salt regime, Eq. ( ^.141 ). The parameters used for e, T and V2 are the 



same as in Figure 11, while i/s = eips/ksT = —0.5 and a = 5A. The bulk concentration 
Cm = 0b is assumed to be much smaller than cm- Adapted from ref. ||112 |. 



Figure 14: The adsorbed amount F as a function of the salt concentration Cgait, Eq. ( p. 18 



for / = 0.01 and 0.25. The solid curves on the right hand side correspond to the scaling 
relations in the high-salt regime, Eq. ( |5.18| ). The horizontal lines on the left hand side 
mark the low salt values, Eq. ( |5.14|) . The dashed lines serve as guides to the eye. The 
parameters used are: e = 80, T = 300 K, V2 = 50A^, a = 5A, yg = —2.0. Adapted from 
ref. [ITT^ . 



Figure 15: Schematic diagram of the different adsorption regimes as function of the 
charge fraction / and the salt concentration Cgait- Three regimes can be distinguished: (i) 
the low-salt regime D <^ (ii) the high-salt regime (HS I) D ^ k,~^ for weak PE's 
f f* = (csaitf2)^'^^; and (iii) the high-salt regime (HS II) D > k,^^ for strong PE's 
/ » /*■ 

Figure 16: Schematic PE brush structure. In a) we show the weak-charge limit where 
the counter-ion cloud has a thickness H larger than the thickness of the brush layer, L. 
In b) we show the opposite case of the strong-charge limit, where all counter- ions are 
contained inside the brush and a single length scale L ^ H exists. 

Figure 17: Scaling diagram for PE brushes on a log- log plot as a function of the graft- 
ing density p and the fraction of charged monomers /. Featured are the Pincus-brush 
regime, where the counter-ion layer thickness is much larger than the brush thickness, the 
osmotic-brush regime, where all counter-ions are inside the brush and the brush height 
is determined by an equilibrium between the counter-ion osmotic pressure and the PE 
stretching energy, and the neutral-brush regime, where charge effects are not important 
and the brush height results from a balance of PE stretching energy and second-virial 
repulsion. The power law exponents of the various lines are denoted by numbers. 
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